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Abstract. In this paper, we present a joint cross-correlation analysis of the Archeops CMB maps at 143 and 217 
GHz and the WMAP CMB maps at 41, 61 and 94 GHz with sky templates of projected galaxy density constructed 
from the 2MASS Extended Source catalog. These templates have been divided in patches sorted in decreasing 
galaxy density with a fixed number of pixels (we considered patches having 32, 64 and 128 pixels) and the cross 
correlation has been performed independently for each of these patches. We find that the densest patch shows a 
strong temperature decrement in the Q, V, W bands of WMAP and in the 143 GHz channel of Archeops, but 
not in the 217 GHz channel. Furthermore, we find that the spectral behavior of the amplitude of this decrement is 
compatible with that expected for the non-relativistic thermal Sunyaev-Zel'dovich effect, and is incompatible (at 
4.5(T level) with the null hypothesis of having only CMB, noise and a dust component (oc u^) in those pixels. We 
find that the first 32-pixel sized patch samples the cores of 11 known massive nearby galaxy clusters. Under the 
assumption that the decrement found in this patch is due entirely to the thermal Sunyaev-Zel'dovich effect, we 
obtain an average comptonization parameter for those sources of i/ = (0.41 ± 0.08) x 10~^ at 13 arcmin angular 
scales. This value is compatible at 1-a with the expectati on, y = 0.49 x 10 ~'^, from a model of the cluster flux 
number counts based on the standard A-CDM cosmology llXue fc Wu 200 J) . The observed value of y is slightly 
diluted when computed for the first patch of 64 and 128 pixels, presumably due to the inclusion of less massive 
clusters, and the dilution factor inferred is also compatible with the quoted model. 
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1. Introduction 

The thermal Sunvaev-Zel'd ovich effect (hereafter tSZ, 
jSunvaev & Zel'dovic h imh) constitutes a unique tool to 
explore the presence of baryons in the Universe. It arises 
as a consequence of the distortion that the black body 
spectrum of the Cosmic Microwave Background (CMB) 
radiation experiences when it encounters a hot electron 
plasma. In this Compton scattering, electrons transfer 
energy to the CMB radiation, generating an excess of 
high energy photons and a deficit in the low energy tail 
of the distribution. This photon reallocation translates 
into a frequency dependent change of the brightness 
temperature of the CMB, which, in the non-relativistic 
limit, has a very simple form, {ftsz{x) = xj tanh[a::/2] —4, 
with X = hv /ksTcMB the adimensional frequency in 
terms of the CMB temperature monopole Tcmb)- 
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The amplitude of this distortion is proportional to 
the electron pressure integrated along the line of sight 
{STtsz/TcMB = ftsz{x) J dr arne ksTe/ {mf,c^), with 
ctt the Thomson cross section and rie, Tf. and mg the 
electron number density, temperature and mass, respec- 
tively); and this makes this effect particularly sensitive 
to collapsed or collapsing structures containing hot 
electrons, such as clusters and superclusters of galaxies 
(see (,Birkinshaw iggQ.,) for a extensive review). 

In addition to the intrinsic energy inhomogeneities 
of CMB radiation generated during inflation, there are 
further temperature anisotropies introduced in the CMB 
during recombination, which are mainly caused by two 
physical processes. These processes are the last Doppler 
kick exerted by electrons via Thomson scattering just 
before recombining, and the subsequent gravitational 
redshift experienced by CMB photons as they climb 
the potential wells generated by the inhomogeneities 



2 



C. Hernandez-Monteagudo et al.: Spectral Indications of tSZ in Archeops and WMAP CMB Data 



in the matter distribu tion (Sachs & Wolfe effect, e.g. 

While all this happens at 
z ~ 1100, a similar scenario can take place at much lower 
redshifts: as the first stars reionise the universe, new free 
electrons are produced which again scatter CMB pho- 
tons, partially blurring primordial anisotropics generated 
during recombination and introducing new ones at much 
larger angular scales. Also, if ^Ia is non-zero, the decay 
of gravitational potentials in linear scales introduces a 
net blueshift in the CMB radiation at late epochs (z < 
2), wh ich is known as the Integrated Sachs Wolfe effect, 
(ISW l|Sachs fc Wolfe 19671) 1. Despite the fact that most 
of these phenomena introduce temperature fluctuations 
of amplitudes larger than the tSZ effect, the particular 
frequency dependence of the latter should enable its 
separation. Whereas th e first generation of CMB ex- 
periments like COBE ijSmoot et al.l99l|) and Tenerife 
ijCutierrez et al. 200Cll aimed to simply detect the largest 
CMB temperature anisotropies in the big angular sca les, 
experiment s like, e.g., Boomerang jMauskonf et al.200n.) , 
ySA te.ubino-Martm et al. 2003* ), Archeops 
l)Benoit et al. 2002|) and WMAP (..Bennett et al. 2003|) 
have already reached the sensitivity and angular resolu- 
tion levels required to probe relatively weak signals like 
the tSZ effect. 

In this work we perform a combined anal- 
ysis of Archeops and WMAP CMB data, 
searching for spectral signatur es of the tSZ 

effect. Previous works llBennett et al. 20031 

Hernand ez-Mon teagudo &: Rubino-Mart m 200 4 
Mvers et al.2004 iHernandez-Monteagudo et al. 2004: 
Afshordi et al. 2004 iFosalaha et al. 200.1: 

Fosalaba fc Gaztanaga 20041) have claimed the de- 
tection of tSZ in WMAP data at different significance 
levels. However, all those studies were exclusively based 
on spatial cross-correlations of large scale structure cata- 
logues with CMB data. In this work, we take advantage 
of the frequency coverage provided by the combination of 
Archeops and WMAP experiments in order to include 
an analysis of the frequency behavior of a signal which 
is spatially correlated with regions hosting large galaxy 
overdensities. 

The sketch of the paper is as follows: in Section [51 we 
summarize the outcome of the Archeops and WMAP 
experiments, whose data products are analyzed as ex- 
plained in Section 121 Section 0] shows our results, which 
are compared with those obtained from WMAP data. 
Their implications are discussed in Section [3 Finally, we 
conclude in Sectional 



2. The Archeops and WMAP data set 

The Archeops (jBenoit et al. 2002|) ^ experiment was de- 
signed to obtain a large sky coverage of CMB temperature 

^ see ,http : //www . archeops . org) 



anisotropies in a single balloon flight at millimeter and 
submillimeter wavelengths. Archeops is a precu rsor to 
the Planck HFI instrument ((Lamarre et al. 20031) . using 
the same optical design and the same technology for the 
detectors, spider-web bolometers, and their cooling 0.1 K 
dilution fridge. The instrument consists of a 1.5 m aper- 
ture diameter telescope and an array of 21 photometric 
pixels operating at 4 frequency bands centered at 143, 217, 
353 and 545 GHz. The two low frequencies are dedicated 
to CMB studies while high frequency bands are sensitive 
to foregrounds, essentially to interstellar dust and atmo- 
spheric emission. Observations are carried out by spinning 
the payload around its vertical axis at 2 rpm. Thus the 
telescope produces circular scans at a fixed elevation of 
~ 41 deg. The data were taken during the Arctic night 
of February 7, 2002 after the instrument was launched by 
CNES from the Esrange base near Kiruna (Sweden) . The 
entire data set covers ^ 30% of the sky in 12 hours of 
night observations. 

For the purpose of this paper, we concentrate in the low 
frequency channels at 143 and 217 GHz. Maps for each of 
the bolometers have been produced from the Archeops 
processed and foreground cleaned timelines, using the 
Mirage optimal map making code llYvon & Mavet 200^ 
as discussed in l|Tristram et al. 20051) . The maps for the 
4 most sensitive bolometers at 143 GHz were combined 
into a single map at 143 GHz and equally the two most 
sensitive bolometer maps were combined at 217 GHz. The 
CMB dipole is the prime calibrator of the instrument. The 
absolute calibratio n error against the dipole as measured 
by COB E/DMR llFixsen et al. 19941) and confirmed by 
WMAP ijBennett et al. 2003|) is estimated to be 4% and 
8% in temperature at 143 GHz and 217 GHz respectively. 
These errors are dominated by systematic effects. The 
noise contribution in the combined maps at 143 and 217 
GHz was computed using Monte Carlo simulations. For 
each bolometer, by using the power spectrum of the noise 
in the time domain data set, we produced fake timelines 
of Archeops noise. These were processed and projected 
into maps following the same procedures used for the 
Archeops data themselves as described before. 



The WMAP satellite mission was designed to mea- 
sure the CMB temperature and polarization anisotropies 
in 5 frequency bands, 23, 33, 41, 61 and 94 GHz with 
a full sky coverage. The satellite was launched in June 
2001 and its first results after the first year of observa- 
tions l|Bennett et al. _2D03^ included the CMB tempera- 
ture and temperature-polarization cross-correlation power 
spectra, as well as full sky temperature maps for each 
of the frequency bands. In this paper we consider only 
data from the high frequency channels, Q [41 GHz], V[61 
GHz] and W [94 GHz], since only for these bands there 
are foreground clean maps available at the LAMBDA site 
: / /IzLinbda . gsf c . nasa . gov/. 



|http 
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3. The Statistical Analysis 

In this Section we outline the correlation analysis per- 
formed on Archeops and WMAP data and the 2MAS S 
Extended Source Catalog (XSC, l|Jarrett et al. 2003tl ) 
on the sky region covered by Archeops. This 
analysis is essentially identical to that applied in 
I Hernandez-Monteagudo et al. 20o4) . paper to which we 
refer for a more detailed description of the statistical 
method. It consists of a pixel-to-pixel comparison of 
Archeops CMB data with a template of the large scale 
structure built from the 2MASS XSC catalog. 

The 2MASS XSC catalog contains approximately 1.6 
million galaxies detected in the infrared filters I, J and K, 
and covers the whole celestial sphere. Those frequencies 
are particularly insensitive to dust absorption, and for 
this reason this catalog can trace the extragalactic struc- 
ture at very low galactic latitudes. The galaxy templates 
built from it take into account the spatial distribution 
of the galaxies and the instrumental bea m of the CMB 
experiment. By using the HEALPix^ (Cor ski et a,1 199<l> 
tessellation of the sky, we built a map with the same 
resolution parameter {Nside — 512) than the one used in 
Archeops and WMAP data. Every pixel was assigned 
a value equal to the number of galaxies present in such 
pixel. The resulting map was then convolved with a 
window function corresponding to the instrumental beam 
of each of the detectors taken into consideration. For 
Archeops, the resulting templates were then weighted 
by their noise levels and co-added per frequency band, in 
such a way that we ended up with two different galaxy 
templates corresponding to the 143 GHz and 217 GHz 
bands. For WMAP, we produced templates for the Q, V 
and W bands. 

In the next step, we sorted the pixels of each tem- 
plate in terms of its amplitude, so that first pixels would 
have higher galaxy densities. These pixels were grouped in 
patches of varying sizes (32, 64 or 128 pixels per patch), 
and again patches were sorted in such a way that first 
patches contained larger projected galaxy densities. Next, 
we analyzed each of these patches separately, by compar- 
ing them to the corresponding patches in the Archeops 
and WMAP CMB maps on a pixel-to-pixel basis. 

As explained in e.g. 

ijHernandez-Monteagudo fc Rubino-Martin 2004|) . it 
is possible to estimate the contribution of a given spatial 
template (M) on a total measured temperature map (T), 
which is the result of the addition of several components: 



J- cmh 



• M + N, 



(1) 



namely CMB (Tcmfc), instrumental noise (N) and some 
signal coming from the extragalactic template (M). The 
contribution of M is parametrized by a, and an optimal 



value of it [optimal in terms of the temperature model 
given in eq.(0), together with its formal error bar, is given 

by 



1 



(2) 



The matrix C is the covariance matrix of T, which must 
contain the correlation matrices of both the CMB and 
noise components. For the small scales we are probing 
here (the pixel has a size of '--^ 7 arcmins), the noise is 
the main contributor to the covariance matrix^. In our 
case, this matrix must be evaluated only for the pixels 
belonging to the patch under analysis. Since our patches 
are relatively small, the inversion of this matrix poses no 
numerical problem. 

Therefore, for every patch a value of a and a a was ob- 
tained. However, uncertainties in the determination of the 
noise amplitude may bias our determination of CTq , and for 
this reason we computed a different estimate of the uncer- 
tainty of a, namely the r.m.s variation of this parameter 
for all available patches, which will be denoted as cr™^ . 
We remark the fact that, according to eq|21 an error in 
the noise normalization does not affect the estimates of 
a, but only those of a a- When comparing a a with u™^ , 
we found that for the 143 GHz channel of Archeops, the 
latter was about a 40% larger than the former, and hence 
we decided to adopt it in order to quote conservative esti- 
mates of statistical significance. For the 217 GHz case no 
such bias was found, and we decided to use again cr™*. 

4. Combined results for Archeops and WMAP 

In figure Q we plot the recovered a's versus the patch 
index. Results are grouped in six different panels: top 
and bottom panels refer to WMAP and Archeops 
experiments, respectively, whereas left, middle and right 
panels display results for patches of 32, 64 and 128 pixels, 
respectively. For Archeops, filled circles and diamonds 
refer to 143 GHz and 217 GHz respectively, whereas 
for WMAP those symbols correspond to the W and 
the Q channels, being the results of the V band given 
by the crosses. Dark and light grey colored bands limit 
the 2-(T confidence levels for filled circles and diamonds, 
respectively, while the moderately dark bands refer to 
the V band (crosses) in the case of WMAP. As explained 
above, for the 143 GHz and 217 GHz channels the 
amplitude of the shaded regions was computed from the 
typical dispersion of the values obtained for a in patches 
where the tSZ contribution is expected to be negligible, 
i.e., in patches with indexes between 40 and 300. We 
have found that for the 143 GHz channel the first patch 
contains an unusual negative a while at 217 GHz seems 
to be compatible with zero. This patch, in the case it 
contains 64 pixels, hosts the central pixels of 20 different 



2 HEALPix's URL site: 
|http : //www. eso . org/science/healpix/ j 



^ Note that for these small scales, the assumption of a pios- 
sonian noise is a good approximation 
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Fig. 1. a's obtained by means of eq.||2J) from Archeops and WMAP data in the fraction of the sky covered by the 
former experiment. In panels (a-c), filled circles, crosses and diamonds refer to the W, V and Q bands respectively. 
In panels (d-f), the 143GHz and 217GHz channel results are given by filled circles and diamonds, respectively. The 
colored bands display the 2-a confidence limit for the a's: dark grey for filled circles, intermediate grey for crosses, 
and light grey for diamonds. Our analysis are focused on the first patch, hosting the highest projected galaxy density. 



AGO clusters of galaxies, COMA among them. Out of 
them, 11 (GOMA again included) are already sampled 
by the 32 densest pixels. Its statistical significance is 
slightly bigger for patches with 64 pixels (> 2.5-(t''™*), 
since it contains the first two very negative patches of 32 
pixels each, (see figure^). In figure^, from the first 
300 patches, very few of them 12) depart from zero 
by an amount similar to that of the first patch; and such 
number is very close to what one expects under Gaussian 
statistics at 2.5-cr level of significance. This peculiar 
behavior of the first (or first two) patch(es) disappears 
at 217 GHz (see diamonds): in no case the diamonds 
corresponding to the first two patches trespass the 
limit. This picture is consistent with part of the signal 
being generated by the tSZ effect, since such compo- 
nent is negative at 143 GHz and becomes zero at 217 GHz. 

In order to interpret the results from WMAP, one 
musts keep in mind that the Q and V bands have remark- 



ably larger beams than the W band: while the Archeops' 
and W band's Point Spread Functions are similar in size 
(« 13 arcmins), the beams of the Q and V bands have an 
average (linear) size of ~ 31 and 21 arcmin respectively. 
This, in terms of tSZ flux, corresponds to factors '--^ 5.7 and 
2.6 smaller in the Q and V band for point-like objects. On 
the other hand, it is clear from figure H that WMAP has 
a much lower noise level when compared to Archeops, 
(approx. a factor of 3.5). As one looks at the top panels 
in figure ^ one finds that for the first patch of both 32 
and 64 pixel size, the W band gives a decrement about 
4-cr away from zero, (which however is not found for the 
second patch of 32 pixel size). This statistical singularity 
of the first patch decreases remarkably in the Q and V 
bands, but still remains close to the 2-a level. However, 
due to the argument on beam dilution on bands Q and 
V with respect to W, we still find the amplitudes given 
by the three bands of WMAP consistent with being (at 
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Fig. 2. tSZ fitted amplitude, A, as a function of patch 
number. We only observe a significant temperature decre- 
ment {A > 0) in patch number 1. 

least partially) generated by unresolved objects causing a 
temperature decrement. 

It also worth to remark that, if the first patch con- 
tains pixels which do not contribute considerably to the 
tSZ decrement, the significance of the overall a obtained 
for that patch will diminish accordingly. This is our ex- 
planation for the decreasing significance of the values of a 
in the first patch for a size of 128 pixels when compared 
to a size of 64 pixels. 

5. Discussion 

5.1. Spectral dependence of the cross-correlation 
coefficients 

As discussed in the previous section, Fig. indicates a 
significant temperature decrement in patch 1 for both the 
Archeops and WMAP data sets which seems to be com- 
patible with the tSZ effect. For a better assessment of this 
result we have compared, via a linear fit, the observed 
spectral dependency of the correlation coefficients a in 
each of the patches to the following model 



a{iy) ^ Ax ftsz{v) + n{v) 



(3) 



where A is global calibration factor which is estimated 
from the fit, ftsz{i^) represents the spectral behavior of 
the tSZ-induced change in brightness temperature and 
n{v) is the instrumental noise in a{v). For this model if 
the signal observed is compatible with the tSZ effect we 
expect a temperature decrement A to be significantly pos- 
itive, otherwise we expect no correlation and therefore A 
to be compatible with zero. We must remark that, due to 
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Fig. 3. The a correlation coefficient in CMB tempera- 
ture units as function of frequency in GHz for patch 1 
(top panel) and patch 201 (bottom panel) of the 32 pixel 
patches. The diamonds are the WMAP data points and 
the stars the Archeops one. In red, we overplot the best- 
fit tSZ computed as discussed in the text. 



the different size of the beams for every channel, we had 
to rescale the a's to a common reference beam size of 13 
arcmin. When doing this, we assumed that the signal was 
generated by unresolved sources, and hence we scaled the 
a's as the ratio of the area of the beam with the reference 
one. 

The main results of the fitting procedure described be- 
fore are shown on figure |21 where we trace the fitted am- 
plitude, A, of the tSZ signal as a function of the patch 
number for patches of 32 pixels. We observe that only for 
patch number 1, where most massive clusters might be, 
there is a significant temperature decrement { A > ). For 
a closer look to the fits, figures 13 01 [S] represent the a coef- 
ficient in CMB temperature units as function of frequency 
in GHz for the 32, 64 and 128 pixel patches respectively. 
The top plot corresponds to patch 1 and the bottom one 
to a patch containing an average value of the projected 
galaxy density (patch 201). In the three figures we ob- 
serve that the correlation coefficients found for the patch 
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Fig. 4. a correlation coefRcient in CMB temperature units 
as function of frequency in GHz for patch 1 (top panel) 
and patch 201 (bottom panel) of the 64 pixel patches. 
The diamonds are the WMAP data points and the stars 
the Archeops one. In red, we overplot the best-fit tSZ 
computed as discussed in the text. 
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Fig. 5. The a correlation coefficient in CMB tempera- 
ture units as function of frequency in GHz for patch 1 
(top panel) and patch 201 (bottom panel) of the 128 pixel 
patches. The diamonds are the WMAP data points and 
the stars the Archeops one. In red, we over plot the best- 
fit tSZ computed as discussed in the text. 



1 are consistent which what we expect for tSZ emission 
with A values of 112± 22 /xK, 77 ± 16 and 54 ± 13 /xK 
respectively. Notice that there is good agreement between 
the model and the data with reduced values of 3/4, 
3.4/4 and 2/4 for patches containing 32, 64 and 128 pix- 
els respectively. For the null hypothesis {A ~ Q /iK), the 
reduced values are 28/5, 26/5 and 20/5. We interpret 
these results as a spectral indication of the measurement 
of a tSZ signal in patch 1. However, for all other patches 
we find A compatible with zero, showing that the data are 
in good agreement with the null hypothesis. For example, 
for patch 201 the A values are 5 ± 19 /^K, —4 ± 13 /iK and 
10 ± 10 with reduced values of 2.2/4, 1.2/4 and 
1.3/4 respectively. We consider these results compatible 
with a no detection of tSZ signal in patch 201. 

The uncertainties for the total amplitude of the tSZ 
signal. A, presented above do not account for systematic 
errors on the Archeops data. These are dominated by 
residuals from atmospheric and Galactic dust emissions 



which in a first order approximation increase as in the 
Archeops frequency range. To account for those contri- 
butions we have recomputed the total amplitude for the 
tSZ signal adding an extra term in the fitted function as 
follows: 



a{v) ^ Ax ftsz{v)+ Asys-xv^ +n{v), (4) 

where Agyg is computed only from the Archeops data. 
For patch 1 we obtain A values of 110 ± 22 /iK, 74 ± 16 
/iK and 53 ± 12 /iK, and Agys values of —95 ± 102 /iK, 
-104 ± 71 /iK and -46± 51 /iK for patches of 32, 64 and 
128 pixels respectively. From these results we conclude 
that the A coefficients are not significantly affected at l-cr 
level by systematic effects present in the Archeops data. 
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5.2. Determination of the mean and integrated 
Compton parameters 

It is interesting to check whether pixels in patch 1 
correspond to potential sources of tSZ or not. For the 
32-pixel-size patches, patch 1 includes COMA, A0576, 
A0671, A0952, A1795, A2061, A2065, A2244, A2245, 
A2249, A2255. Since these are massive and relatively 
nearby galaxy clusters, it is reasonable to expect some 
signature of the tSZ effect. Using the quoted values for 
the A parameter we infer the following estimate for 
the average Compton parameter in all those sources: 
y = (0.41 ± 0.08) X 10-". In the case of 64-pixel-size 
patches we have roughly the same clusters and we expect 
therefore the signal to be diluted, y = (0.28±0.06) x lO"'^. 
Finally for 128 pixel-size patches, patch 1 includes the fol- 
lowing 27 clusters : COMA, A0077, A0104, A0272, A0376, 
A0407, A0576, A0671, A0952, A1035, A1185, A1235, 
A1377, A1767, A1795, A1800, A2034, A2061, A2065, 
A2069, A2142, A2151, A2199, A2244, A2245, A2249, 
A2255. We deduce for them y = (0.20 ± 0.05) x 10"^ 
Since we have scaled the a's to the beam-size of the 143 
GHz channel of Archeops, these quoted values of y are 
associated to a (linear) angular scale on the sky of ~ 13 
arcmins. For the above results we ha ve assumed the CM B 
temperature to be Tcmb = 2.725K ijMather et al. 1999(l . 

We now try to relate the observed average tSZ decre- 
ment {yobs) to the high end of the SZ number counts. We 
make 2 hypotheses: 1) that the Archeops beam encloses 
most of the integrated tSZ effect in clusters and 2) that 
the 2MASS survey is a perfect tracer of the tSZ effect. We 
will come back to these hypotheses at the end. 

The correlation analysis, that is presented in the previ- 
ous section, can be recast in stating that the Npi^ brightest 
clusters of galaxies have an average integrated Y = J dily 
parameter eq ual to Y = {yobs)^beam where flbeam = 
27r(FWHM/\/81og2)2 is the Archeops beam sohd an- 
gle if FWHM = 13 arcmin. The integrated Compton 
parameter can be directly related to the SZ flux num- 
ber counts whi ch has been the issue of a great number 
of studies (e.g. ("Aghanim et a l. 1997tlBartlett et al. 19941: 
[Barbosa et al. 1996; Bartelm arm^tal. 200ll)) from w hich 
we select llXue fc Wu 200ll) a^^Benson et al. 2002|) . 

Following f Xiie fc Wu 200lh there are three possible 
models for which the number counts of clusters over the 
whole sphere can be parametrized in terms of y as follows, 



N{>Y)^No{Y/Yo)-\ 



(5) 



where we fiducially consider Yq = lO^^arcmin^ 
For Model 1 (Ml), deduced from the cosmological 
ACDM matter power spectrum, A'o — 29 and 7 — 1.5. For 
Model 2 (M2), a non-evolving X-ray luminosity function is 
used to correct the counts and gives a larger Nq ~ 635 with 
the same exponent. Finally Model 3 (M3), an evolving X- 
ray luminosity function is used instead, giving iVo = 350 
with a flatter exponent 7 = 1.2. The conversion from flux 



in Jansky units to the Compton frequency-independent 
quantity Y is obtained via 



-rr — J^CMB Jtsz[J^)- 



(6) 



^arcmin? is 



For example, the fiducial value Yq — 10 
equivalent to 0.75 and 0.91 Jy at 90 and 143 GHz (resp.). 
From the above formulas, we deduce the lower limit 



Yrnin in the number counts for the first N^i 
brightest clusters as 



Yjy^^j^f^^ Yf) 



-1/7 



' Nofsky ' 



32,64,128 



(7) 



where fsky ^ 0.20 is the effective sky fraction used in the 
Archeops- WMAP tSZ cross analysis. The average Y value 
of these clusters is then 



{Y)^Yo 



7 



-1/7 



1-1'Nof 



sky 



(8) 



We deduce the average Compton parameter for those 
clusters (y) as 



(y) 



(Y) 



27r(FWHM/V81og2)2 



(9) 



For the previous three models, expected values for (y) 
are between 0.5 x 10"'* and 6 x 10^'*. The value expected 
from basic principles (Ml), (y) = 0.49 x lO^"', is quite 
close to the observed value {y)obs — 0.41 ± 0.08 x IQ-"*. 
The observed value should however be corrected up- wards. 
Non-linearities in the relation between the 2-MASS den- 
sity field and the Y parameter introduce an efficiency 
which is difficult to estimate, although we note that as 
discussed above most of the brightest pixels in the density 
map constructed from the 2-MASS survey are associated 
to well-known clusters of galaxies. If we considered that 
only identified clusters (11) are present in the first patch of 
32 pixels, the above expected values should be multiplied 
by a factor of two. Extension of the bright clusters beyond 
the fiducial beam of 13 arcmin (like Coma) may produce 
a differential effect with the observing beam at different 
frequencies, as well as a bias in the number counts. 

We also note that the dependence of (y) with the num- 
ber of pixels taken in the analysis (32, 64 or 128) follows 
quite well the Model 1 predicti on as shown in Table ^ 
We have also cross- checked the dXue fc Wu 200lh model- 
ing with that from Ib enson et al. 2002^ . This alternative 
model produces an intermediate value between models Ml 
and M2 and is marginally compatible with the observa- 
tions. 

6. Conclusion 

In this paper we have presented a joint cross-correlation 
analysis of the Archeops and WMAP data sets with a 
template of galaxy density constructed from the 2MASS 
XSC galaxy catalog. We have first divided the 2MASS 
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sky density map in patches of equal number of pixels 
and sorted these patches in terms of decreasing projected 
galaxy density. For each of these patches we have per- 
formed a cross-correlation analysis with the Archeops 
data at 143 and 217 GHz and with the WMAP data 
for the Q, V and W bands. For patches containing the 
densest 32, 64 and 128 pixels (patch 1), the correlation 
test pointed to a prominent temperature decrement in 
WMAP's Q, V and W bands and in the 143 GHz band 
of Archeops, but not at 217 GHz, as would be expected 
for tSZ-induced temperature fluctuations. All the other 
patches failed to show a similar behavior. 

To corroborate these results, for each of the patches 
we have compared the cross-correlation coefficients to a 
model of the tSZ frequency pattern, in which we fix the 
spectral behavior and fit for a global amplitude parameter 
A. For the first patch and the three sizes considered (32, 
64 and 128 pixels), we obtain non zero A values at more 
than 4.5-iT level with good agreement between the model 
and the data, and negligible contribution from systematic 
residuals in the Archeops data. For all other pixels hav- 
ing smaller projected galaxy density, we fail to detect any 
signature of tSZ effect. 

From these results, we conclude that there is clear in- 
dication of tSZ effect for patch 1 in the Archeops and 
WMAP data sets. This is not surprising, since patch 1 
contains pixels centered in massive and relatively nearby 
galaxy clusters. Assuming the signal observed is tSZ, we 
infer, for the 32-pixels case, an average value for the comp- 
tonization parameter oi y = (0.41 ± 0.08) x lO""* in all 
those clusters at an angular scale of ~ 13 arcmins. This 
value is compatible at l-cr level with the expectations, 
y = 0.49 X lO"'^ (cf. Table HJ, from a model of the clus- 
ter flux number counts based on the standard A-CDM 
model. Ml, assuming the measured y is due to the con- 
tribution from the 32 brightest clusters. For 64 and 128- 
pixeles patches the tSZ signal is diluted, probably due to 
the contribution of relatively not so massive clusters. Note 
that the dilution observed is also compatible with the one 
expected from the model Ml (cf. Table 0. 
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Npix 


CMB Data 


Ml 


M2 


M3 


B 


32 


(0.4 ±0.08) 


0.49 


3.91 


6.24 


2.02 


64 


(0.27 ±0.06) 


0.31 


2.46 


3.50 


1.27 


128 


(0.19 ±0.04) 


0.19 


1.55 


1.96 


0.80 



Table 1. Expected average Compton parameter, (y) /lO in patch 1, as a function o f the number of pixel s in the patch 
for the four number counts models discussed in the text, Ml, M2 and M3 l|Xue k Wu 2001,) and B (jBenson et al. 2002f) 



